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We study the primordial magnetic field generated by stochastic currents produced by scalar 
charged particles created at the beginning of the radiation dominated epoch. We find that for 
the mass range 10~ 6 GeV < m < 10 2 GeV, a field of sufficient intensity to seed different mecha- 
nisms of galactic magnetic field generation, while still consistent with observational and theoretical 
constraints, is created coherently over a galactic scale. 
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Q\ ' I. INTRODUCTION 
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At present there exists huge observational evidence about the presence of magnetic fields throughout the Universe: 
our own galaxy is endowed with a homogeneous magnetic field B ~ 3 x 10~ 6 Gauss and similar field intensity is 



< 



detected in high redshift galaxies Q and in damped Lyman alpha clouds 

The origin of these fields is still unclear. Research is performed mainly along the idea of a cosmological mechanism: 
a seed primordial field that would be further amplified by protogalactic collapse and differential rotation or a nonlinear 
dynamo ||. Several mechanisms have been proposed to explain the origin of the seed field. It has been suggested 
that a primordial field may be produced during the inflationary period if conformal invariance is broken |jj , |5j . In 
string-inspired models, the coupling between the electromagnetic field and the dilaton breaks conformal invariance 
and may produce the seed field ||. Gauge invariant couplings between the electromagnetic field and the space-time 
^vq curvature also break conformal invariance but produce, in general, an uninterestingly small seed field 0. Other 
£■ — ■ mechanisms are based on a first order cosmological phase transition || and on the existence of topological defects || . 

In this paper we propose a new mechanism for primordial field generation in the early Universe, based on stochastic 
currents generated by particle creation of scalar charged species |0L We assume the presence during Inflation of a 



charged, minimally coupled, scalar field in its invariant vacuum state [by. When the transition to a radiation dominated 
Universe takes place, quantum creation of charged particles occurs. We assume that the field mass is smaller than the 
Oh' vacuum energy density during Inflation, H, and therefore can consider the transition from that period to Radiation 
q , dominance as instantaneous. 

5^ ■ As the number of positive charged species is the same as the number of negative charged ones (there is no physical 
reason why it should not be so), the mean electric current is zero. Nevertheless quantum fluctuations around the 
mean give rise to a non- vanishing current. We compute the rms amplitude of these fluctuations and use it as the 
source term in the equation for the magnetic field. We must stress the fact that the field must be a scalar minimally 

• <-h , coupled one: it is straightforward to check that with a massive conformally coupled scalar field very few particles are 
created and therefore a very weak magnetic field is created. For spinorial fields we show in the Appendix that, due 
to the conformal invariance of massless fermions, the number of particles created is very small and consequently the 

- - ■ magnetic field produced is extremely weak. 



II. CALCULATION OF THE MAGNETIC FIELD 



As the process of magnetic field generation that we are studying takes place after inflation, there is no loss of 
generality if we work in a spatially flat Universe. In conformal time, dr — dt/a(t), we have g^ v = a 2 {r)rj^, r\^ v being 
the minkowskian metric tensor. The canonical scalar field, which we assume to be massive and minimally coupled, is 
written as ip = 4>/a. 

Part of this process takes place before the electroweak transition. Then part of the created photons are a combination 
of of the isospin and hypercharge bosons, where the coefficients are respectively the sin and cos of the Weinberg angle 
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9 W H , |T^] . We will consider the magnetic field generated by only the hypercharge sector and in this sense the figures 
obtained constitute a lower bound to the effective intensity of the field. The amplitude of the electromagnetic field 
would be smaller than the pure U(l) boson by a factor cosO W: but recalling that the electroweak coupling constant is 
g = qj cos9 w , we obtain the same amplitude for the created electromagnetic field. The magnetic field is then defined 
from the spatial components of the field tensor, Bi = {\/2)eijkFjk, where F^ u = d^A v — d^A^, being the vector 
potential in the Lorentz gauge A^ = 0. 

The equation for B reads: 



9 2 „2 / x 9' 



B = Vxj (1) 



where a (r) is a time-dependent conductivity of the plasma and j represents the spacelike components of the electric 
current which for the scalar field is = iq ((f)*d^(f) — cfyd^*) + 2q 2 A fJj (p*(f>. 

The basic point in our analysis is that while the expectation value of the current vanishes, there are quantum and 
statistical fluctuations which build up a nonvanishing root mean square value. The source in these equations can 
be phenomenologically substituted by a stochastic, Gaussian current p3| . Our goal is to estimate the magnetic field 
produced by this stochastic source. Under the assumption that the source is Gaussian, its statistics are fixed once 
the two point source-source correlation function is given. Since the radiated field will be weak in any case, as a first 
approximation we can compute this correlation function in the absence of a macroscopic electromagnetic field. 

As = 0, the only contribution to the electric current will be due to the quantum fluctuations of the fields. In 

order to study if the quantum fluctuations of the massive fields produce a non zero electric current we calculate the 
two point correlation function 

N iV [(r, e({(vx J) >,r), (V x j) v (r',f)}) (2) 

The scalar field can be written in terms of two real fields according to </>(t, r) — {4>i(t, f) + i<j>2{T, r}} /V2 so that 
the spacelike part of the current j M , at vanishing e.m. field reads j — q |</>2V(/>i — </>iV(/>2 j and 

V x j = 2g(v0 2 ) x (V^i) (3) 

If we consider that 4>\ and </>2 commute, we can write 

N W [(r,r), (tV)] = {d? f G+df lt G+ + g>?.,G"c$,G-} (4) 

where we have introduced positive and negative frequency propagators G + = (</>i(r, r)</>i(r', r 1 )) and G~ = 
(</>i(r' ,r')(f>i(T,r)) . Writing the propagators in terms of their Fourier components we get 

JV«' l(r, r), (r\ ?)] = V J x k'\ (k x k'\, e^'^-^G^ r')G+ (r, r') + (r, r < — > r', f) (5) 

Expanding the fields as (j) = (2tt)~ 3 ^ 2 j d 3 n K (r)e IK '" + h.c. the Fourier transform of the propagators is given by 

Gt(T,T>)=cj) K (T)^(T') = G-(T>,T). 

Our next task is the evaluation of the scalar field modes. In the absence of electromagnetic fields, the scalar field 
satisfies the Klein - Gordon equation 



9 , ..2 , 22/ n a(r) 



-(- k -\- m z a' i {T) 



dr 2 a (t) 

We assume instantaneous reheating at r = 0, 



Mr) = (6) 



Inflation 



H(t -t) 

<>lr) = { (7) 
■jpp ( 1 + ^- J Radiation 

and normalize the scale factor by taking Htq = 1. In terms of the dimensionless variables yo = Ht and k — k/H the 
modes in the inflationary epoch read 
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bl(yo) = 4V^VoH^ [k(l - y )] , v = %\\\-\% 



For the radiation dominated era we write 



9H 2 



(8) 



where /P is the WKB solution 



% KB (yo) = a k f k +p k f k * 



(9) 



e -iCl k (y ) 



where Wfc(yo) = \J k 2 + jji-a(yo) 2 , ^fe(yo) = J V ° dy' uj k (y' ). a k and (3 k are the Bogolubov coefficients connecting the 
modes (0) with the WKB basis fM) at yo = 0. In the limit k <C m/H <C 1 they are given by 



Pk — —otk — —i 



r(3/2) ttH 1 



fc 3/2 



(10) 



The non zero value of (3 k reflects particle creation from the gravitational field, rather than from the decay of the 
inflaton field. Indeed the occupation number for long wavelength modes diverges as k~ 3 , much in excess of the k~ x 
Rayleigh - Jeans tail of the thermal spectrum produced by the reheating process. This excess noise results in an 
enhancement of the magnetic field over and above the equilibrium fluctuations. When the given scale enters the 
horizon, particle - antiparticle annihilation becomes efficient, and the extra noise disappears. It is worth mentioning 
that the logarithmic divergence in the total number of particles created up to fc max can be removed with a suitable 
infrared cut-off, say the mode corresponding to the horizon ~ 10~ 26 . The energy density associated to these particles, 
p = mH 3 J m ** d 3 k \Pk\ 2 , after a cutoff is imposed is ~ 10~ 16 smaller than the one of the CMBR. 

To give a quantitative estimate of the magnetic field amplitude, let us come back to the evaluation of the two point 
correlation function (jsj) , that can be written as Nw = N^, + N^, + Nf it . The first term is the noise that would be 
present in the absence of particle creation during reheating and does not concern us here. The other terms are the 
contribution due to the particles created at the beginning of the radiation era and are given by 



Nl, [(r,r),(r',r-')]^4g 2 ^ 



dk dk' 



kxk')(kx k') i e<^')-(y-y)G lk (yo,y' a )SG lkl (y Q ,y' ) 



(11) 



Nl, [(T,r),(r'y)] -> q 2 H*J ^ (k x (k x k') 6^')^') SG lk (y , y' )6G lk , (y , y' ) (12) 

where 

G° lk = G+° + G- k ° = C - ^MBM (13) 

v^fc(yoVfc(2/o) 

with Q K (yo, y' ) = ^«(yo) - M K (y' ), 

SG lk (y , y' ) = 2a k (3if k (y )fM) + ^ k fk* (Vo) f° K * (y'o) + 2 \Pk\ 2 G lk (Vo,y'o) (14) 

The energy density of the magnetic field can be calculated from the two point correlation function generated by the 
stochastic current: 

(B(x)B(x>)} =H± J ^^^G ret (y ,y;x ,x)G ret (y' ,g;xo,x') x (15) 



x {Nh [(y , y) , (y' , jf)] + N 2 [(y , y) , (y' , y>)}} 

Where G ret are the causal propagators for the Maxwell field, obtained from the homogeneous solutions to Eq. (Q). 
By transforming Fourier the causal propagators, the spatial integrals can be readily performed by virtue of the simple 
spatial dependence of the propagators. We have 
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(B(x)B(x')) = H 4 / dy dy' x 



dk dk' 



kxk 1 



x {4<?G° lk ( y(h y' Q )6Gi k ,(y ,y' ) 



(16) 



The energy density of the magnetic field coherent over a given dimensionless spatial scale A is given by E\ — 
-^y? J dx J dx" (B(x)B(x' )) which amounts to insert the window function W kk '(X) = V^ 1 J v ^dx e l \ k + k ')- x . The 
energy density reads pb = a~ A (B 2 ), where 



(Bl) = H 4 



kxk' 



GJk t +k'\(yo,xo)G^ t +k , l (y' ,x ) 



\k+k'\\ 



dyody' x J |^W fe V(A) 
x {4q 2 G lk {yo,yo)6G lk '(yo,y' ) + q 2 6G 1 k{yo,yo)SG lk '(yo,y' )} 



(17) 



k + k' 



A < 1 and W fc fe'( A ) ~ for 



k + k' 



A > 1, we can take as 



As the window function satisfies W kk >(\) « 1 for 

the upper limit in the momentum integrals, fc m ax — j and W k k' = 1 in the interval (0, fc max ). All the cosmological 
interesting scales are such that k <C 1 : For example for a galaxy we have that its physical scale today is (if it had not 
collapsed) A G ~ IMpc ~ lO^GeV- 1 = X c T rh /T today ~ 10 28 A C (T rh ~ 10 15 GeV is the temperature of the Universe 
at the end of reheating and Ttoday — 10~ 13 GeV is the present temperature of the microwave background). The 
commoving wavelength is A c ~ lO^GeV -1 , taking H = lO n GeV the dimensionless scale A = HX C ~ 10 21 and the 
corresponding dimensionless momentum fc ~ 10~ 21 . 

The retarded propagators are constructed with the homogeneous solutions of Eq. (|l|) . If we take the conductivity 
(cfr. Refs. g|) a ~ Te~ 2 = T rh e- 2 (1 + yo^ 1 = a (1 + yo)" 1 , 
a = a /H = 



-T 



rh 



I pi 



He 2 y / m pl /H, 



larger than the horizon) 



! \J rripi /H being the dimensionless conductivity, the homogeneous equation (Q) reads (for scales 



jB k + 



c 



dyl (1 + yo ) dy 



B k =0 



(18) 



where we have used B k (yo) = J d 3 ye lk y B(y, yo). The causal propagator for this equation is given by 

ff — l" 



/~<ret 
^\k+k>\ 



i + yo 
l 



CO 



1 - 



1 



Vo 



1+Xq 



(19) 



If xq 3> yo the final expression for the propagator is GT^f.,, ~ (1 + J/o)/co- 

In spite of the fact that the scalar field is not in thermal equilibrium with the background radiation, the electro- 
magnetic interaction with it causes a correction to the value of the field mass, proportional to the temperature of 



the bath, i.e. we have m 2 



eT 2 



e ^r/i/(l + Vo) ■ We will consider this correction in the frequencies 



w fe(j/o) but not in the Bogolubov coefficients because we assume instantaneous particle creation. In fact, when the 
transition from Inflation to radiation dominance occurs, large numbers of of particles are created very quickly by 
parametric resonance. Thermalization is a process that takes place afterwards and during a longer period of time [ pT[ . 
As for the ratio m/H, if we consider, for example, the Higgs mass (to ~ l0 2 GeV), we have m/H ~ 10 -9 . We can 
therefore neglect the terms k 2 in front of those m 2 /H 2 . When we perform the products G® k (yo, y' )SGi k ' (yo, y' ) and 
&G\ k (yoi y'o)SGik' (yotUo) m Eq. (0) we w ^ nave terms where the exponentials cancel and terms where this does not 
occur. The latter oscillate and hence can be neglected because they will give a smaller contribution than the first ones 
when the time integrals are performed. These integrals are performed between and T, where T is the dimensionless 
time when a scale reenters the horizon; after this time, the current quickly relaxes to its equilibrium value through 
particle - antiparticle annihilation. From the expression of the Bogolubov coefficients, Eqs. ( |l0| ) we can see that the 
contribution from the term quartic in (3 k overwhelms the quadratic one. The leading contribution to ( |l7|) reads 



H 4 k 4 — 




ml (l+j/o)_ 
if 2 



eT 2 



H 2 



q 2 = e 2 = 1/137 and recalling T 



21 , for the electroweak mass scale, too — 10 2 GeV, ao — e 2 10 4 , T r h ~ 10 
7r/ (2/e max ) we first note that, by simply replacing these figures in Eq. (p 



(20) 

15 GeV, 
), the 
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term proportional to m 2 , (1 + yo) overwhelms the one proportional to T 2 h in the square root by a factor ~ 10 20 . 
Therefore we can take 

(Bl) ~ \f (—) 6 H*k^^ = \? (-) Vfc^ x ^ (21) 
3 \m J erg H z 3 \m J at 



(Bl) ~ 10 24 GeU 4 -> r = — ~ 1(T 36 (22) 

Pfocfe 

where p6 C fc = T 4 ^ is the energy density of the CMBR. These values correspond to a commoving field of strength 
B\ ~ 10 32 Gauss. The physical field = a~ 2 B\, such that = B 2 h , gives a present value of 

B ph ~ 10" 24 Gauss (23) 

This value satisfies the constraints imposed by the anisotropy in the microwave background Jlq ] and Big Bang 
nucleosynthesis |l6| and (marginally) suffices to seed the galactic dynamo (|J^] , [0. Stronger fields are obtained 
by diminishing the value of the field mass. We can estimate the lower bound to the mass by demanding that today 
B p h ~ 10~ 9 Gauss, according to the mentioned constraints. We obtain m > 10~ 6 GeV. Considering the value of the 
mass used throughout the paper as the upper bound we find the range 10~ 6 GeV < m < 10 2 GeV. Fields outside of 
this range will contradict observational and/or theoretical constraints. 



III. CONCLUSIONS 



To conclude, we have presented a new mechanism for the generation of a primordial magnetic field, based on 
the breaking of conformal invariance. In previous works conformal invariance was broken in order to amplify 

quantum fluctuations of B. The amplification was in general very small. Here scalar, massive charged particles 
minimally coupled to gravity, and coupled to the electromagnetic field produce stochastic fluctuations in the source 
of Maxwell eqs. which in turn generate an astrophysically relevant primordial field. For conformally coupled massive 
scalar field, the Bogolubov coefficient scale with the momenta as k" 1 / 2 which means that the number of created 
particles is very small and consequently the magnetic field generated will be too weak to be astrophysically interesting 
(~ 10~ 80 times smaller than in the minimally coupled case, for the physical parameters used in this paper). The 
fact that spinor fields are conformal invariant if massless, manifests itself in that the proposed mechanism does not 
generate strong enough magnetic fields, as is shown in the Appendix. Finally we must say that the value of the field 
quoted in Eq. (^3|) is to be considered as a rough estimate; a more precise evaluation will require a more detailed 
consideration of the whole electroweak gauge theory and of the non-equilibrium evolution of current and field. 



IV. APPENDIX 



Here we show that the magnetic field generated by an electric current originated by the creation of fcrmionic 
particles is indeed negligible. The Dirac equation for spinors in a FRW Universe reads (in conformal time) ]il||ll| 

[i-y^-ma(yo)]x(yo,y) = (24) 

where x(yo,y) = « 3 ^ 2 (2/o)'0(j/O; V)- Ho and y are the dimensionless time and space coordinates defined in the paper 
and in this Appendix m is to be understood as a dimensionless mass, i.e. m — > m/H. The electric current reads 

j" = ex(yo,2/Vx(2/o,y) 

where xiVOiV) = xHyo,y)l° , i- e - the Dirac adjoint. The positive and negative energy spinors, Uk s and Vk s read 

^ V vyu; Wj fk(yoV h - y (25) 



Uks 



Vks = 



C s 
C s 



r k (yo)e- rk -y (26) 
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where fk(yo) is a solution to the equation 



d . ,2 , 22/ \ ■ 9a(y ) 
oVq dy 



A = 



(27) 



The spinors are normalized according to the (time independent) product J" s dT,xj°X 
The noise kernel due to this field reads (after a rather long calculation) 

AT t , 2 f ,r f ,Jk{y'o)f P {y'o)fk(yo)fp(yo) 
Nu(yo,y Q )=^ jdkjdp NkNpNkNp 



X7°X? S being a spacelike hypersurface. 



where 



[M>o,^)+A<S(2/o,^) + 



fc +p 



^S(yo,yo) = -^ 



ma(y ) 



■ Am 

/fc(j/o) 



ma(y ) 



/ \ /fc(yo) 

™a(y ) - i ; 

4(yo) 



■Ami 

Wo). 

/;(!»)' 



, x ../pU/o; 



(28) 



(29) 



k. (k + pj (k + pj .p 



k 2 p 2 
ma(y ) - i 



ma(y' Q ) + i 



Jk(y'o) 



fk(y' Q ) 



ma(y' ) 



Ml/o) 

f P (y'o) 



(30) 



A* 0/q) 
fk(yo) 



ma(y ) - i 



j;{yo) 



f;(y ) 



Nk, p are the normalization factors which we calculate for Inflation. The solution to the field equation (g7j) for the 
inflationary period that corresponds to positive frequency for yo — > — oo read 



f k (yo) = (1 - yo) 1 ' 2 HP [k (1 - yo)] ; v = \ + im 



and therefore 



N k = N k 



1 



For the radiation dominated period and for k/m <C 1 we write the mode functions as ppf 

/>oo 

h = a fe e mi(1+ ^ )2/2 + p k e lm ^+y^ 2 / 2 / e- s2 ' 2 ds 

Jz{ya) 



(31) 



(32) 



(33) 



where z(yo) = (1 + i)^/m(l + yo) and a k and j3k are the Bogolubov coefficients obtained by matching the modes and 
its time derivatives at yo = 0. In the limit k/m -C 1 they read 



A=-i r B +i -Mfy /2+ '™* 2 



(34) 
(35) 



We see that in this case there is no divergence in the number of quanta created and therefore the expected electric 
current will be very small. We replace the modes (|3^) in the equations ( p9~l) and (^) and with the obtained expressions 
we evaluate 



(B 2 )=H± I ^^G^\yo,xo)G^(y' 07 xo)N u (y ,y' )Wl p (X) 

J (27T) (27T) 



(36) 



G 



with the same considerations that were made for the scalar field. The main contribution comes from the first two 
terms and they are respectively 

f ,s Q p 2jm(l+T max ) 2 

(B 2 f ] - 2 _ 2[ 2]2 TT. , C (37) 

m o 1 1 ' m \ cosn [itm] 

(B 2 f ] -q 2 ^— — % = J&„ (38) 

to 2 (Tq [j + m 2 J cosh [nm] 

For example, we have that for a mass m = 10~ n (i.e. physical mass ~ f GeV in units of H = 10 11 GeV), fc max and 

a the same as the ones used in the paper, (B 2 ) {1) ~ 1CT 160 GW 4 and (B 2 ) {2) ~ lCT^GeV 4 , values completely 
negligible in comparison with the one quoted in Eq. (p2). 
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